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ABSTRACT 


'  Jji 

Tlag_oi|^od--empioyed“  by  •The-Dow-Chemical-Company'  for  calculating 
specific  Impulses  of  propellant  mixtures  . is-cle8ori-bed-j'--'Al-th©ugh 
method-l^general-,  »l^he  capacity  of  the  computer  (Burroughs 
220)  and  the  availability  of  thermodynamic  data  place  practical 
restrictions  on  what  systems  can  be  calculated.  The  calculations 
center  around  the  method  of  White,  Johnson,  and  Dantzig  for 
determining  the  equilibrium  compositions  of  complex  mixtures. 

By  developing  and  applying  matrix  equations  which  describe  both 
the  theraodynamic  data  and  the  physical  equations,  the  program- 
mlng  effort  has  been  reduced  .signlfj^q^anjtly.  Required  input  data 
to  the  program  are  the  composition  of  the  reaction  mixture,  the 
chamber  and  exit  pressures,  and  a  rough  estimate  of  the  chamber 
temperature.  All  other  pertinent  data  are  retrieved  from  magnetic 
tape  by  the  computer  program.  The-program-is"'amoma1ri{r*IiT'TRat''" 
^large  number  of  problems  may  be  solved  without  manual  inter¬ 
vention.  In  addition  to  the  specific  impulse,  the  program 
calculates  the  characteristic  velocity,  the  chamber  and  exit 
temperatures,  compositions,  enthalpies  and  entropies. 
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I.  INTRODUCTION 


The  theoretical  specific  impulse  of  a  propellant  system  has  been 
used  extensively  as  an  indicator  to  evaluate  the  maocimum  per¬ 
formance  which  could  be  expected  from  the  system.  The  calculation 
of  the  specific  impulse  allows  rapid  and  relatively  inexpensive 
screening  of  the  many  candidate  systems.  Earlier  work,  previously 
reported^^®,  at  this  laboratory  by  R.  R.  V.  V/iederkehr  resulted 
in  a  general  method  programmed  for  a  Burroughs  220  Computer. 

Recently,  the  capabilities  of  the  Dow  Burroughs  220  were  enlarged; 
therefore,  expansion  of  the  computational  program  was  undertaken. 
The  basic  mathematical  model  remained  unchanged,  but  the  setup, 
input,  and  output  sections  were  extensively  revised. 

Although  the  basic  method  has  been  discussed  in  the  previously 
cited  reports,  for  the  sake  of  unity  and  clarity  an  exhaustive 
treatment  of  the  theory  and  modus  operandi  of  the  program  will 
be  included  in  this  report. 

Specific  impulse  has  been  defined  as  the  velocity  of  the  products 
at  the  exhaust  divided  by  the  gravitational  constant.  If  the 
change  in  enthalpy  of  the  chamber  products  as  they  are  expanded 
through  a  nozzle  is  equated  to  the  kinetic  energy  of  the  particles, 
the  specific  impulse  may  be  defined  as: 


v/here,  is  the  enthalpy  of  [J-  grams  of  products  in 
the  rocket  chamber, 

Hg  is  the  enthalpy  of  ki  grams  of  products  in 
the  nozzle  exit,  and 

c  is  a  constant. 


In  order  to  calculate  the  specific  Impulse  one  must  merely 
determine  values  of  Hq  and  Hg.  Assuming  no  heat  losses  to  the 
surroundings,  a  heat  balance  in  which  the  enthalpy  of  the  react¬ 
ants  is  converted  to  the  enthalpy  of  the  products  may  be  written: 


r' 

I 

j=i 


=  H 


(2) 
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v/here,  is  the  number  of  moles  of  the  jth  reactant 
In  the  reaction  mixture, 

is  the  heat  of  formation  of  the  jth 

reactant  at  298 “K.,  and 

p  is  the  number  of  reactants. 

assumed.  Prom  this  Jxhaust  Pee" 

be  easily  calculated  as  a  funct?on  r>r°^  enthalpy,  Hg,  can 

sary  Information  for  tL  tenperature.  All  the  neces- 

equatlon  (1)  Ttn°n  ‘P®  PPP'^lflo  l>"pulse  by 

u“de?groSo1?  cl:jn\e'”''Sge"egeni°or'«  ,  ’’''“‘'““‘p 

ascertained  by  theoretical  mean?*^  Tt  change  cannot  be 

eSSanslon^J^M^S'"^”^^  equilibrium  Is  malntalned^duSlng 

ca!  prlS?e^'’S^ch"'a“''mereJv'’tiff?  p'PP'-pPp 

equlllbrlum  composuIon  S  a  aljL  ^  PhP»'leel 

desired  are  fixed  bv  the  enthfimr  temperature.  The  temperatures 
chamber  solution  and  by  the 

exhaust  solution  as  illustrated  in  Table  I 

Table  I 


Chamber:  H(T^,  X^,  Pc )  =  He 
Exhaust:  S(Te,  Xg,  )  =  Sc 


He  and  Igp 


tempemLjf  and^cLpoJJtlof  freach^'of  tJe  twfv^^hl 

rx^\^on^\f^bTL^n:rrf?rlal\^^df^- 

left  side  of  the  equation  ^  P^°°®dure  and  then  the 

temperature  and- composition,  if  equallt-^°ls^ob? 

equality  is  not  obtained  an  in:proved1eSlSure°irL^^imel  "" 
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The  Magnitude  of  the  Improvement  is  dependent  on  the  magnitude 
and  direction  of  the  inequality.  The  equilibrium  temperature  Is 
again  calculated  and  the  process  repeated  until  equality  is 
attained. 

It  Is  evident  that  the  solution  to  the  problem  revolves  around  a 
method  for  determining  the  compositions  and  temperatures.  This 
report  will  consider  the  following: 

(l)  Method  for  calculating  composition  (Section  II). 

(ll)  Method  for  calculating  temperatures  (Section  III). 

(ill)  Method  for  generating  thermodynamic  functions 
(Section  IV) . 

(iv)  Description  of  the  operation,  input,  output,  and 
capabilities  of  the  computer  program. 


II.  CALCULATION  OF  THE  EQUILIBRIUM  COMPOSITION 
OF  A  rWLTICOMPONENT  SYSTEM^ 


Prom  the  first  and  second  lav;s  of  thermodynamics  it  can  be  shown 
that  a  system  at  constant  temperature  and  pressure  achieves 

energy  of  the  system  is  a  minimum. 
Equilibrium  is  achieved  v/hen  the  composition  is  such  that  the 
free  energy  is  minimized.  This  principle  is  used  to  calculate 
the  equilibrium  compositions  of  multicomponent  systems. 

The  method  used  for  minimizing  the  free  energy  of  a  system  is 
essentially  that  of  V/hlte,  Johnson,  and  Dantzlg®  as  extended  by 
Thomson*  to  Include  condensed  species*. 

The  method  performs  the  following  operation:  given  any  composi¬ 
tion  Y  =  (yi,  ya,  ya  . .  .y^)  compatible  v;lth  the  mass  balance  con¬ 
straints,  an  Improved  composition  X  =  (xi,  Xg,  X3  ...x^)  is 
calculated  which  is  closer  to  the  equilibrium  value.  If  X  Is 
not  sufficiently  close  to  the  equilibrium  value,  it  becomes  the 
new  Y,  which  can  be  Improved  further.  The  procedure  is  repeated 
until  convergence  is  achieved. 


*This  section  is  taken  almost  in  its  entirety  from  Report  Nr. 
AR-lS-60,  by  R.  R,  V.  Wlederkehr  (Reference  (l)). 
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A.  Composition  Convergence  Method 


The  free  energy  of  the  system  with  a  composition  of  X  is*: 
nt 

^  f/PA 

+  In  P  +  In 


Pr(X) 


=  t  [(i) 

i=l  ‘- 


X  +  y  ^ 

X  i  A  RT 

i=n„+l 


"g 

where,  x  =  ^  . 

i 


tig+ne  ^ 

V  Pi 

2»  BT  *1 

(5) 

l=ng+l 

(^) 

The  summation  in  (4)  extends  over  only  the  gaseous  components. 

The  first  summation  on  the  right  of  (5)  is  the  free  energy  of  the 
gaseous  species;  the  second  term  is  the  free  energy  of  the  con¬ 
densed  species.  The  entropy  of  mixing  of  the  condensed  phases  is 
assumed  negligible,  and  the  gases  are  assumed  to  be  ideal. 

The  behavior  of  the  free  energy  in  the  neighborhood  of  Y  may  be 
approximated  by  a  Taylor's  expansion,  about  Y: 


ng+n^ 


Q(X) 


=  Pr(y)  +  X 


JJL 


i=l 


(x^-yj^) 


X=Y 


"g+"c  '^g+^c 

I  X  ^ 


•p 


(x^-y^)  (Xj^-yj^)  (5) 


X=Y 


which  reduces  to: 

n. 


Q(X)  =  F^(X)  + 


X  [( 

'PA  Yi' 

^  )  +  In  P  +  In 

1=1 

y  J 

(x^-yi) 


g  c 


x'  ®  I 

i=l  1 


i=ng+l 


(6) 


*P^(X)  is  actually  free  energy  divided  by  RT. 
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It  Is  desired  to  minimize  the  quadratic  approximation  to  the  free 
energy  Q(X)  subject  to  the  element  balance  restraints  given  by: 


Z  =  ‘■j 

1=1 


(7) 


vjhere  a^j  are  formula  numbers  indicating  the  number  of  atoms  of 
element  j  in  a  molecule  of  species  i,  and  bj  is  the  total  number 
of  atomic  vjeights  of  element  J  originally  present  In  the  mixture; 
m  Is  the  number  of  elements  in  the  system. 

The  constraints  and  Q(X)  may  be  incorporated  into  a  single 
equation: 


0(X)  =  Q(X)  *  (v  +  bjJ 

Where  the  tTj's  are  Lagrange  multipliers. 

G. 


(8) 


To  minimize  G(X) , 


is  set  equal  to  zero.  It  can  be  shovm 


via  (6)  and  (8)  that  for  the  gaseous  components: 


_J_G 


\B.T  J 


+  In  P  2  In  — i- 


+ 

1  X 

Yx  "  y 

m 


- 

J=i 


1=1,2, 


.n 


g 


^ij  =  0 


(9) 


but  for  the  condensed  components: 


i  G 
i  x^ 


P, 

RT 


m 


jaij  °  ^“"g+r  \+2’  •••"g‘^"c 


(10) 


There  are  ng  equations  of  type  (9),  n^  of  type  ( 10) ,  m  constrain¬ 
ing  equations  (7),  and  equation  (4).  At  this  point,  therefore, 
there  are  n^  +  m  +  1  equations  and  as  many  unknowns .  But  by 
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jS/Pf  s  and  substituting  these  x^'s 
equations  (7)  and  (4),  the  fiumber  of  unIaio\flis  Is  reduced  tc 
Uq  +  m  +  1.  Solvlrus  (9)  for  yields: 


into 


4  =  -yi^i  + 


(w-a  V,) . 


1=1,  2,  ...ng  (11) 


v/here 


1  RT/ 


+  In  P  +  In 


Substitution  of  (ll)  into  (4)  and  (7)  produces: 

m  "g  "g 

I  X  ’"j  =  X  ^1^1  ( 

j-1  Li=i  J  “i 

V  Tr®  1  1  - 

I  L  \  ^  X  I  *  X  ^ij"i  = 

K  i  Li-1  J  ji-2  J 


e> 

'’j  *X 


The  problem  of  ralnlralzlng  the  free  energy  subject  to  mass  balance 
constraints  reduces  to  solving  equations  (l4),  (lo),  and  (13) 
simultaneously. 

shown  that  these  equations  can  be  rewritten  in  matrix 
notation  as  follows : 


\d  V  y-  ir  =  B  +  A'.  Y, 


A^tt  =  y 

C  J’C 


C  C  g  g  y 


gCi'g 


''■g  V  -  5 


g  ^  g 
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where  the  matrices  are  defined  as  follows 


Ag 


ail  ai2> 

asi 


®im 

I 

I 


I 


Ac 


ang+1,1*  •  •  ang+i,m" 

^  • 


'  •  J 

Ung+ncji*  •  'ang+ncni 


ng  X  m 


Hq  X  m 


Pi' 

yi 

h 

ya 

Rt 

• 

1 

1 

• 

Png 

Yg  » 

1 

1 

1 

- 

^gD  =• 


yx  0  0  0 

I 

0  ya  0  I 

0  0  ya  ! 

!  0 

_6  ....  OVn^ 


Xg 


Xi 

fa 

xng 


B  = 


bm 

The  superscript 


bi 

ba 


Xc  = 


ng+i 


Lfng+ncJ 


designates  the  transposition  of  a  matrix 


Equations  (l5)j  (l6)  and  (17)  may  be  combined  to  give  a^slngle 
matrix  equation: 


m  nc 


■ 

/  ■  -  A  ■  ■■  ■ 

r* 

m. 

^gD  ^ 

A'c  I  Ag'Yg 
1 

TT 

nc. 

• 

Ac 

0  *  0 

1 

“Xc 

Y'g  Ag 

0  •  0 
•  . 

/ 

/  -  — '> 

B  +  A'g  YgD  3g 


(18) 


This  equation  represents  a  system  of  m  +  n^,  +  1  linear  equations 
with  as  many  unknowns.  The  brackets  and  symbols  Indicate  the 
dimensions  of  the  matrices. 


The  unknowns  are:  ir  =  (Vu  iTg,  ...Tr|p)f 

^c  =  (^ng+i>  ^ng+2  *  • ‘^ng+ng)  > 


Solving. Equation  (l8)  for  the  unknowns  yields: 


•  • 

TT 

A'gYgDAg;  A'elVXg 

1  1 

Xc 

__  A,  ;  0  ;  0 

"x'g'Ag”:  "oT’o' 

-r 


®  ^'g  ^gD^g 


y 


g  Jg 


(19) 


The  square  matrix  to  be  Inverted  may  be  calculated  from  Y  and 
the  formula  niombers  a. ..  The  vector  on  the  right  may  be  calcu¬ 
lated  from  Y,  a.^,  ana'^the  3-p.  andj-^,  vectors.  It  will  be  shown 
In  Section  IV-A^now  the  3-  and  ^  vectors  (partial  free  energy 
vectors)  may  be  simply  calculatea.  Therefore,  given  Y,  the 
formula  numbers,  and  the  partial  free  energy  vectors,  one  may 
calculate:  the  number  of  moles  of  each  condensed  species  at 
equilibrium,  Xq,  the  Lagranglan  multipliers  tt,  and  x/y.  The 
number  of  moles  of  each  gaseous  species  at  equlllbrlvim  may  be 
calculated  from  Equation  (ll)  .  In  matrix  notation  Equation  (ll) 
Is 

r  ^  — *  ~i 

(20) 


^  -  ’'ei’Kf )  -  ''g”] 


where  both  Xg  and  [Cf)  ■  AgT^are  (ng  x  1^ ^vectors,  and 

•g  X  ng)  diagonal  matrix. 


YgD  Is  a  (ng 


In 


ii-  A  V 

(20)  (I) 


represents  a 


(ng  X  l)  column  matrix  whose  elements  are  all  x. 

y 


-8- 


At  a  givsn  temperature  the  thermal  data  and  the  formula  numbers 
Hence,  Equations  il9)  and  (20)  provide  a  scheme  for 
obtaining  an  improved  estimate  of  the  equilibrium  composition, 
fpom  an  original  estimate  Y.  The  vector  Y  may  be  regarded 
as  an  input  to  our  scheme  which  produces  an  output  X.  The 
resulting  X  becomes  the  nev;  Y  and  the  procedure  is  repeated 
until  X  and  Y  balance. 

It  is  not  practical  to  continue  the  procedure  until  an  exact 
balance  of  X  and  Y  is  obtained,  since  this  would  take  an 
Inordinate  amount  of  time.  Therefore,  the  following  criterion 
is  used; 

If  for  all  yi  [xi  -  yi)^^^  <0.0005  (21) 

Then  X  is  assumed  to  be  in  balance  with  Y. 

B.  Convergence  Difficulties 

Certain  difficulties  may  arise  during  the  process  of  balancing 
X  and  Y,  i.e.,  calculated  values  of  xi  may  be  negative,  the 
square  matrix  in  equation  (9)  may  fail  to  invert,  or  X  may  con- 
tlnuously  oscillate  about  Y  and  never  quite  converge.  There 
follows  a  discussion  of  how  each  of  these  difficulties  may  be 
circumvented. 

1.  Negative  Values  of  xi 

It  sometimes  happens  that  values  of  Xi  calculated  by  the  above 
method  are  negative.  Physically  this  has  no  significance  since 
the  number  of  moles  of  any  species  xi  must  be  zero  or  positive. 
To  rectify  this  condition  when  It  occurs,  the  following  pro¬ 
cedure  Is  follov/ed; 

A  =  X  -  Y  (22) 

which  represents  the  recommended  step  needed  to  Improve  the 
compositions.  If  one  or  more  of  the  Xj['s  are  negative  the 
recommended  step  Is  too  large.  Therefore,  It  Is  reasonable  to 
merely  shorten  this  recommended  step  by  a  factor  A  so  as  to 
make  each  of  the  xi's  positive.  The  new  values  of  xi  are  given 


Y 

new 


Y  +  AA 


(23) 


for  each  x^  which  Is  negative  let  a  Ai  be  calculated  according 
to  the  equation. 


>^1 


-Yl 

“ST 


(24) 
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Substitution  of  this  Xi  Into  Equation  (25)  causes  ^ew 

zero.  If  a  A  larger  than  Xj.  v/ere  used,  xi  new  v/ould  remain 
negative.  Therefore,  of  all  the  calculated,  the  smallest 
must  be  chosen.  This  value  of  X,  viz. .  vrlll  cause  Xs  new 
to  be  zero,  and  all  other  xi's  to  be  positive. 


a.  Gaseous  Species  -  For  gaseous  species,  zero  Is  an 
objectionable  value  to  use  for  Xs  since,  according  to  equation 
(20),  this  estimate  can  never  be  Improved  but  v/111  remain 
zero.  Instead  of  letting  xs  go  to  zero,  the  following  scheme 
Is  used:  Arbitrarily  values  are  given  to  Xg  depending  on  the 
size  of  the  corresponding  ys,  that  Is, 


if 

Ys 

>  1  X  10“^^, 

Xg  =  1  X 

10-11 

(25a) 

If 

1  X 

10"^^  >  yg  >  1  X  lO"^^, 

Xg  =  1  X 

10"  ^1 

(25b) 

if 

Ys 

<  1  X  10"^^, 

Xg  =  0 

(25c) 

If  ys  Is  small  enough,  Xs  does  become  zero.  Since  this  v/ould 
preclude  the  possibility  of  any  further  modification  or  Improve¬ 
ment  to  this  specie,  when  convergence  has  been  attained  all 
zero  values  are  arbitrarily  set  at  a  value  of  1  x  10"°  and 
the  Iteration  repeated  until  X  and  Y  balance  again.  This 
device  Is  a  radical  departure  from  the  previously  reported 
method^-  and  has  proved  to  be  more  efficient.  Every  time  a 
new  temperature  Is  used  the  device. 

If  XI  =  0,  yj[  -=  1  X  10  (26) 

Is  also  used. 


b.  Condensed  Species  -  For  condensed  species  a  zero  value 
for  Xs  Is  not  objectionable.  Consequently,  using  As  alone  Is 
sufficient. 

Another  difficulty  may  arise,  however.  If  xi  for  a  condensed 
species  Is  zero  (which  means  y±  for  the  next  Iteration  Is  also 
zero)  and  the  following  x±  Is  negative,  an  Impasse  Is  reached. 

Clearly  Ai  calculated  from  Equation  (24)  Is  zero  so  that  Xnew  = 
Yoid  [see  Equation  (25)]*  A  A  of  zero  prevents  X  from  Improving. 
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III.  CALCULATION  OF  THE  CHAMBER  AND  EXIT  TEMPERATURES 


Tivo  temperatures,  the  chamber  temperature  and  the  exhaust 
temperature,  must  be  determined  in  the  course  of  the  solution 
of  the  problem.  The  desired  temperature  is  obtained  when  the 
enthalpy  of  entropy  of  the  system  equals  the  chamber  enthalpy 
or  entropy  respectively  (See  Table  I).  Since  the  methods  used 
for  the  tv/o  temperatures  are  similar,  only  the  case  of  the 
chamber  temperature  will  be  considered. 

Starting  with  an  assumed  temperature,  T,  the  enthalpy  of  the 
system  is  calculated  and  compared  v;ith  the  chamber  enthalpy. 

If  the  calculated  enthalpy  is  low  (high),  then  the  next  tempera¬ 
ture  is  increased  (decreased).  This  new  "improved  temperature," 
Tn+i  is  calculated  according  to  the  following  linear  interpola¬ 
tion  or  extrapolation  formula: 


Tn+1  =  Tn  +  MUc  -  Hn)  (27) 

where  Hn  is  given  by  Equation  (55)  and  He  is  given  by  Equation 
(2) .  To  calculate  the  second  temperature  from  the  assumed 
first  temperature  a  value  for  of  0.01  +K./cal.  is  used.  For 
succeeding  temperature  is  given  by: 


Tn"Tn-i 

Hn-Hn-1 


(28) 


Equation  (27) is  used  for  extrapolation  to  new  temperature  until 
the  chamber  temperature  has  been  straddled.  Once  the  chamber 
temperature  has  been  straddled.  Equation  (28)  is  used  for 
Interpolation. 

Equation  (27)  has  certain  shortcomings.  When  the  temperatures 
Tn  and  Tn_i  straddle  both  the  chamber  temperature  and  a  phase 
transition  temperature.  Equation  (27)  produces  a  very  poor 
Improved  temperature."  This  Is  a  result  of  the  very  non¬ 
linear  behavior  of  H  with  T  In  the  neighborhood  of  a  phase 
transition;  the  enthalpy  abruptly  changes  across  the  transition 
(See  Figure  l) .  When  this  difficulty  arises  It  can  be  eliminated 
by  the  following  method  which  Is  based  on  Figure  1. 


Enthalpy 

or 

Entropy 


Fig.  1  -  Schematic  of  Non-Linear  Behavior 
of  Enthalpy  v;ith  Temperature 

Suppose  that  the  tv/o  temperatures  which  straddle  both  the  chamber 
temperature  and  a  phase  transition  temperature  yield  the  points 
A  and  B  on  Figure  1.  By  calculating  the  enthalpy  of  the  system 
at  a  temperature  slightly  (O.S^K. )  above  and  below  the  phase 
transition  temperature,  the  points  C  and  D  are  obtained.  If 
the  chamber  enthalpy  is  betv/een  the  enthalpies  of  points  A  and 
C  or  B  and  D,  Equation  (27)  is  applied  using  the  points  A  and 
C  or  B  and  D,  respectively.  Hov/ever,  if  the  chamber  enthalpy 
lies  betv/een  the  enthalpies  of  points  C  and  D,  then  the  chamber 
temperature  coincides  v;lth  the  phase  transition  temperature  and 
the  compound  undergoing  the  phase  transition  must  be  partly  in 
one  phase  and  partly  in  the  other. 

In  general,  the  composition  of  the  system  varies  upon  traversing 
a  phase  transition.  The  phase  transition  composition  (chamber 
composition)  may  be  calculated  on  the  assumption  that  the  com¬ 
position  varies  linearly  with  the  enthalpy  across  the  transition. 
From  the  enthalpies  and  compositions  corresponding  to  C  and  D 
and  the  chamber  enthalpy, the  chamber  composition  can  be  cal¬ 
culated. 

To  this  point  the  discussion  has  been  concerned  with  choosing 
Improved  temperatures.  When  an  improved  temperature  T  is 
obtained  such  that : 

H(T,  X)  -  He  <  100  cal.  (29) 

then  that  temperature  T  is  defined  as  the  chamber  temperature 
Tc-  In  the  above  equation  H(T,  X)  is  given  by  (55)  and  Hq 
is  given  by  (2).  The  convergence  criterion  for  the  exhaust 
solution  is ; 

If  S(T,  X,  Pg)  -  Sc  <  0.02  e.u.  (50) 
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then  the  temperature  is  the  exit  temperature.  In  this  equation  Sc 
Is  evaluated  according  to  (65)  using  the  chamber  conditions  and 
S(T,  X,  Pe)  In  calculated  by  (65)  at  the  assumed  exit  conditions. 

IV.  THE  GENERATION  OP  THERMODYNAMIC  PROPERTIES 
BY  THE  USE  OF  MATRICES 


A.  The  Matrix  Representation  of  the  Enthalpy,  Entropy,  and 
^ree  Energy^  - 

To  calculate  the  chamber  conditions  It  Is  necessary  to  be  able 
to  evaluate  the  enthalpy  as  a  function  of  temperature,  composi¬ 
tion,  and  pressure.  To  calculate  the  exit  conditions  It  Is 
necessary  to  be  able  to  evaluate  the  entropy  as  a  function  of 
temperature,  composition,  and  pressure.  To  calculate  the 
equilibrium  composition  It  Is  necessary  to  know  the  partial 
free  energies  as  a  function  of  temperature,  pressure,  and  composi¬ 
tion.  This  section  shov/s  how  these  thermodynamic  functions 
can  be  concisely  represented  and  simply  generated  by  the  use 
of  matrices. 


1.  Single  Species  System 

The  enthalpy  of  a  single  species  (l.e.,  a  compound  In  a  given 
phase)  at  a  given  pressure  as  a  function  of  temperature  can 
be  evaluated  If  the  heat  capacity  as  a  function  of  temperature 
and  the  enthalpy  at  some  base  temperature.  Tv,  are  known: 

H(T)  =  HCT’d)  +  /  Cp  (T)  dT  (31) 


The  molar  entropy  of  a  species  as  a  function  of  temperature  at 
a  given  pressure  can  also  be  evaluated  If  the  heat  capacity  as 
a  function  of  temperature  and  the  entropy  at  some  base  tenpera- 
ture ,  Tb ,  are  known ; 


s(t)  =  S 


Cp(t) 

T 


dT 


(52) 


The  molar  free  energy  of  a  compound  can  be  evaluated  according 
to  the  equation; 


P(T)  =  H(T)  -  TS  (T)  (33) 

Therefore,  given  the  heat  capacity  as  a  function  of  temperature, 
and  the  enthalpy  and  entropy  at  some  base  temperature,  Tb,  one 
can  evaluate  the  enthalpy,  entropy,  and  free  energy  as  functions 
of  temperature  at  a  given  pressure. 


*Thls  section  Is  also  taken  In  Its  entirety  from  R.  R.  V.  Wlederkehr, 
Report  No.  AR-ls-60. 
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Haler  and  Kelly*  have  successfully  described  the  heat 
capacity  as  a  function  of  temperature  using  the  equation: 

Cp  -  a  +  bT  +  ^  (34) 

where  a,  b,  and  c  are  constants  determined  from  experi¬ 
mental  data. 

Thompson*  has  arrived  at  a  significantly  Improved 
expression  for  the  heat  capacity,  viz.. 

Cp  »  a  +  bT  +  cT®  +  dT®  +  eT“®  (35) 

More  generally  the  heat  capacity  can  be  represented  by 
a  function  which  Is  linear  with  respect  to  the  parameters, 
vlz.^ 

N 

Cp  =■  E  Ij  fj(T)  (36) 

Where  the  Ij’s  are  constants*  corresponding  to  the  a,  b, 
c’s  and  the‘^fi(T)‘s  are  pure  temperature  functions. 
Equations  (34)  and  (35)  are  specific  cases  of  Equation 


If  (36)  la 

substituted 

Into  (31) 

and  1 

[  32)  then 

H(T)  = 

H(Tb) 

N 

*  I 

fj{T) 

dT 

(57) 

Tb 

and 

N 

fl(T) 

S(T)  = 

S(Tb) 

+  z 

dT 

m 

J=1 

''Tb 

X 

If  the  Integrals 

in  (37)  and  (38)  are 

defined  by: 

Ij(T)  = 

dT 

(39) 

Jj(T)  = 

/ 

■fjW 

T 

dT 

(40) 

*At  a  given  pressure  the  Ij's  are  constant. 
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H(T) 


f  N 

=  H(Tb)  -  IjIj(Tb)  +  IjIj(T)  (41) 


and 


N  N 

S(T)  =  S(Tb)  -  ^  IjJj(Tb)  +  ^ 


(42) 


Now  if  the  following  definitions  are 

Introduced 

N 

^N+l  -  H(Tb)-  } 

ljlj( Tb) 

(a) 

In+1  =  1 

(b) 

(^3) 

Jn+1  =  0 

(c) 

•Hi 

1n+2  =  S(Tb)  -  ^ 
J=1 

Ij  Jj(  Tb) 

(a) 

In+2  =  0 

(b): 

(44) 

'^Ni-2  =  1 

(c)J 

then  (4l)  becomes 

N-f2 

H(T)  =  Ijlj(T) 

(^5) 

and  (42)  becomes 

N-Fa" 

3(T)  =  IjJj(T) 

(46) 

Substitution  Of  (45)  and  (46)  Into  (53)  yields  an 
expression  for  the  free  energy.  O' -Lexus  an 
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P(T)  =  ^  IjKj(T) 


where 

Kj(T)  5  Ij(T)  -  T  Jj(T)  (48) 

The  equations  for  Hj  Sj  and  F  apply  to  a  single  compound 
In  a  single  phase^  l.e.,  a  single  species.  The  li's 
may  be  Interpreted  as  elements  of  a  thermodynamic ''property 
vector,  1,  which  la  Independent  of  temperature.  The 
Ij,  Jj,  Kj  may  be  Interpreted  as  elemen^  enthalpy, 
entropy  and  free  energy  thermal  vectors  I,  j,  K,  which 
depend  only  on  temperature  and  are  Indepen^nt  of  the 
compound.  The  Inner  product  of  1  with  X,  J  and  if  yield 
the  enthalpy,  entropy  and  free  energy,  respectively. 

Before  the  above  Equations  (45)^^  (^7)  .can  be  of 

any  practical  use,  the  vectors  I,  j,  K  and  T  must  be 
evaluated.  All  of  these  vectors  depend  on  the  linear 
form  chosen  to  represent  the  heat  capacity  data 
(Equation  36). 

In  Appendix  I,  these  vectors  are  calculated  using  the 
linear  form  proposed  by  Thomson  (Equation  35)  and  thermo> 
dynamic  property  vectors  are  tabulated  for  all  of  the 
high  temperature  combustion  products  for  which  data  are 
readily  available. 

Multicomponent  Systems; 

For  the  1th  species  In  a  multicomponent  system  In  which 
mixing  effects  are  negligible,  the  enthalpy,  entropy  and 
free  energy  are: 

N+2 

Hi  ■=  Z,  lljlj  ('tS) 


N+2 

I 


lljjj 


Pi  =  t  l^jKj 
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The  total  enthalpy  of  the  system  Is: 
n 

H  .  (52) 

where  Xj^  Is  the  number  of  moles  of  the  1th  species  present 
In  the  system  and  n  Is  the  number  of  species  In  the  system* 
Similar  expressions  hold  for  the  total  entropy  and  total 
free  energy  of  the  system.  Substitution  of  (49)  Into  (52) 
yields: 


n 

N+2 

H 

"  iS 

^1 

I 

j=i 

or  In 

terms 

of  matrices 

• 

• 

H 

=  [xi. 

X2, 

...  Xjj 

] 

111  ll2 

■ 

• » •  1 1 >  N+2 

lai 

t 

1 

1 

1 

y 

i 

1 

ft 

• • •ln>N+2 

(53) 


(54) 


I 

I 


L^N+2 


or,  more  simply: 

H  =  X'  L  I  (55) 

where  X  ,  L  and  I  are  the  matrices  defined  by  (54)  and 
(55).  Similarly,  matrix  expressions  for  the  total  entropy 
and  total  free  energy  of  the  system  may  be  written: 

s  =  X'  L  J  (56) 

(57) 


P  =  X 


L  J 
L  K 


The  matrix  L  Is  seen  to  be  an  array  of  the  thermodynamic 
property  vectors  of  the  species  present  In  the  system. 
Henceforth,  the  matrix  L  shall  be  referred  to  as  the 
thermodynamic  property  matrix. 


Equations  (55)x  (56)  and  (57)  suffer  from  the  fact  that 
mixing  effects  have  been  neglected.  Mixing  can  be 
accounted  for  by  making  the  following  simplifying 
assumptions: 

1.  All  gaseous  species  behave  Ideally. 

2.  Condensed  species  have  zero  enthalpies  and 
entropies  of  mixing.. 
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For  a  mixture  of  Ideal  gases,  the  enthalpy  of  mixing  Is  • 

zero.  The  entropy  of  mixing  of  the  1th  species,  however. 
Is: 


ASi  =  -  R(Pi 


where 


(Pi  =  In  XI  +  In  £ 

X 


(58) 

(59) 


Consequently,  the  corrected  entropy  and  free  energy  per 
mole  of  the  1th  gaseous  species  are: 


Si  =  Si  -  R(Pi 
Pi  =  Fl  +  RT(Pi 


(60) 

(61) 


respectively.  These  corrected  values  are  also  partial 
molar  quantities. 

Equations  (60)  and  (6l)  are  applicable  to  condensed 
species  as  well  as  the  gaseous  species  If  It  Is  under¬ 
stood  that 

<^1  =  ^  (62) 
an  condensed  species.  Therefore,  Equations  (49),  (60) 
and  (61)  may  be  regarded  as  being  general  expressions  for 
the  partial  molar  enthalpy,  entropy  and  free  energy  of 

species  of  a  multicomponent  system.  Substitution 
of  (60)  Into  (50)  and  (6l)  Into  (5l)>  respectlvelj&  pro- 

UUO  6  S • 


iJ 


N 


Si  =  ^  lijJj  -  R(Pj^ 
J  ^ 


and 


N 

^1  =  L  lljKj  +  RTiPi 


(65) 


(64) 


The  results  of  this  section  can  now  be  applied  to  the 
overall  specific  Impulse  calculations. 

To  determine  the  chamber  and  exit  temperatures  It  la 
necessary  to  be  able  to  evaluate  the  enthalpy  and  entropy 
as  functions  of  temperature  and  composition  (the  pressure 
Is  given).  Equation  (55)  provides  a  very  simple  means  of 
evaluating  the  enthalpy. 

H(T,  X)  =  X  LI  (55) 


O  i 
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If  Equation  (56)  Is  corrected  for  mixing  according 
(63),  the  following  matrix  equation  for  evaluating 
the  entropy  Is  obtained: 

S(T,  P,  X)  =  X'[LJ  -  R<P] 


to 

(65) 


In  Equation  (18),  which  Is  used  to  calculate  an  improved 
composition  convergence,  the  vectors  3g  and  3© 
required.  Both  of  these  vectors  can  be  generated  at 
once  by  the  use  of  Equation  (64)  divided  by  RT.  In 
matrix  notation  this  equation  becomes; 


(66) 


The  subscripts  g  and  c  refer  to  the  gaseous  and  condensed 
species,  respectively  and  the  dimensions  of  the  matrices 
are  designated  in  the  equations  by  brackets.  The  vector 
(Pq  Is  zero. 

Equations  (55),  (65)  and  (66)  are  the  main  results  of 
this  section.  These  equations  require  that 'the  matrices 
L,  I,  J,  K,  (P  and  X  are  given.'  Elements  of  the  matrix 
L  are  independent  of  temperature,  pressure*  and  compo¬ 
sition.  For  a  given  problem,  therefore^j^  L  Is  .Invariant . 
The  temperature  determines  the  vectors  I,  T,  K  while  the 
pressure  and  composition  determine  the  vector  (P.  The 
vector  X,  of  course.  Is  the  composition.  Appendix  I  Is 
concerned  with  evaluating  the  vectors  1,  I,  J,  and  K. 

B.  Arrangement  of  the  Product  Species  and  the  Construction  of 
the  ^hermodynaimlc  Property  Matrix,  L~ 

The  matrix  L,  which  Is  composed  of  the  1  vectors,  depends  on 
how  the  T  vectors  are  arranged.  Since  each  1  corresponds  to 
one  of  the  product  species,  the  matrix  L  depends  on  how  the 
product  species  are  arranged. 


In  the  method  for  calculating  the  equilibrium  composition 
described  In  Section  II,  a  strong  distinction  was  made  between 


*The  elements  of  L  are  actually  functions  of  pressure;  but  at 
high  temperature  the  variations  of  the  l^j  with  pressure  are  I 
probably  small. 
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gaseous  species  and  condensed  species.  One  type  of  equation 
was  used  to  describe  the  gaseous  species;  another  was  used  to 
describe  the  condensed  species. 

Because  of  this  distinction,  care  must  be  exercised  in  keeping 
track  of  in  v/hat  state  each  compound  finds  itself.  At  a  tempera¬ 
ture,  T,  all  compounds  whose  boiling  points  are  less  than  or 
equal  to  T  are  gaseous,  v;hereas,  all  compounds  whose  boiling 
points  are  greater  than  T  are  condensed.  Consequently,  as  the 
temperature  of  a  system  varies,  the  number  of  condensed  species 
may  vary. 

A  general  means  of  treating  this  problem  is  to  consider  all 
species,  i.e.,  all  phases  of  all  compounds,  as  being  present. 

Each  species  would  be  either  gaseous  or  condensed,  and  the 
species  could  be  grouped  accordingly.  Furthermore,  these  two 
groups  v/ould  remain  Invariant  during  the  course  of  the  solution, 
regardless  of  the  temperature  and  pressure.  The  physical  absence 
of  a  given  species  at  a  given  temperature  and  pressure  should  be 
a  natural  result  of  the  minimization  of  the  free  energy,  viz. , 
the  given  species  should  have  an  equilibrium  composition  of 
zero  moles. 

This  general  method,  unfortunately,  has  certain  disadvantages. 

The  number  of  condensed  species  v/ould  always  be  unnecessarily 
large  so  that  the  difficulties  described  in  Sections  II,  C  1 
and  2,  v/ould  be  accentuated.  To  overcome  these  difficulties,  it 
v/as  shov/n  that  certain  condensed  species  must  be  dropped. 
Consequently,  it  seems  best  to  consider  only  those  condensed 
species  which  are  physically  present  under  a  given  set  of  con¬ 
ditions,  and  reject  all  those  which  are  physically  absent.  A 
second  disadvantage  is  that  including  all  species  in  the 
problem,  whether  they  be  physically  present  or  not,  unnecessarily 
Increases  the  computer  storage  and  time  requirements. 

The  following  scheme  was  employed  to  distinguish  between  the 
various  phases.  Since  all  the  possible  gaseous  species  are 
considered,  including  the  gaseous  phase  of  the  condensable 
species,  the  phase  rule  only  allows  one  temperature  at  which 
both  the  liquid  and  solid  phase  of  any  one  specie  will  exist. 

This  temperature  is  the  melting  point,  and  for  the  purposes  of 
this  discussion  the  normal  melting  point  will  be  used.  At  any 
given .temperature  there  will  be  present  either  the  solid  or  the 
liquid  phase  of  a  condensable  species  and  the  gaseous  phase. 

The  method  proceeds  along  the  following  lines:  All  the  compounds 
are  arranged  into  two  groups,  the  gaseous  and  those  which  are 
condensed. 
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It  is  convenient,  therefore,  to  form  tv,'o  matrices  at  the  beginn¬ 
ing  of  a  problem:  the  Lj.  matrix  and  the  L^*  matrix.  The 
matrix  is  composed  of  rows  corresponding  to  the  r  always  gaseous 
compounds;  it  is  invariant  during  the  course  of  the  problem.  The 
Lu*  matrix  contains  tv.'o  rows  for  each  condensable  conpound, 
one  row  for  each  of  the  liquid  and  solid  phases.  At  a  given 
temperature  one  rov/  for  each  conpound  is  taken  from  the  Ln* 
matrix  and  used  to  form  the  Lq  matrix  which  characterizes  the 
thermal  properties  of  the  condensable  compounds  at  that  tempera¬ 
ture.  The  normal  melting  point  of  each  compound  determines 
which  row  of  Ly*  is  to  be  used  in  Lq. 

The  matrix  L  v/hich  characterizes  the  thermodynamic  properties  of 
the  entire  system  is  composed  of  the  tv/o  submatrices  Lr  and  Lq. 


V.  THE  COMPUTER  PROGRAM  FOR  CALCULATING 
THE  SPECIFIC  IMPULSE 

A.  Flow  Diagram  for  Calculating  the  Specific  Impulse 

Figure  2  is  a  flow  diagram  which  describes  both  the  overall 
procedure  for  calculating  the  specific  impulse  and  the  computer 
program  based  on  this  procedure.  The  overall  procedure  may 
be  subdivided  into  three  parts: 

(i)  Data  handling  and  preliminary  calculations 
(ii)  Chamber  calculation 
(iil)  Equilibrium  exhaust  calculations 

Each  row  of  Figure  2  corresponds  to  one  of  these  parts.  Before 
the  methods  for  calculating  temperatures  and  compositions  can 
be  employed,  certain  data  must  be  available.  These  data  must 
be  selected,  retrieved,  and  calculated,  and  as  Is  schematically 
illustrated  by  the  first  row  of  Figure  2 .  These  computer 
performed  operations  are  discussed  below. 

1.  Read- in  Data 


Input  data  consist  of  the  composition  mixture,  the  chamber  and 
exit  pressures,  a  rough  estimate  of  the  chamber  temperature,  and 
any  necessary  Instructions  needed  to  modify  the  number  of  pro¬ 
duct  species  considered.  (See  Appendix  II) 
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2  -  Composition  Convergence  Method 


2.  \fliat  Elements  are  Present? 


A  number  of  predetermined  reactants  are  stored  along  with  their 
code  number,  molecular  formulas,  heats  of  formation  at  298®K. , 
molecular  weights,  names,  and  Identifying  words  consisting  of 
the  appropriate  atomic  numbers.  This  Identifying  number  con¬ 
tains  the  atomic  numbers  of  all  the  elements  of  the  compound  in 
ascending  numerical  order. 

3.  Construct  Reactant  Matrices 

For  each  reactant,  data  are  retrieved  from  the  magnetic  tape, 
the  Ap  and  C  matrixes  are  constructed,  and  a  list  of  each  of 
the  following  is  compiled:  molecular  weights,  heats  of  forma¬ 
tion,  and  print  words  for  the  molecular  formulas  and  names. 

4.  What  Products  are  Possible? 

A  number  of  predetermined  products  are  also  stored  with  their 
code  numbers,  formulas,  thermodynamic  data,  molecular  weights, 
melting  points,  and  identifying  numbers  such  as  those  described 
for  the  reactants.  Products  v;hlch  contain  only  those  elements 
present  in  the  fuel  mixture  are  possible  species.  This  Is 
accomplished  by  comparing  the  atomic  numbers  of  the  elements 
present,  found  in  part  1  of  this  discussion,  v;ith  the  atomic 
numbers  In  the  Identifying  word.  If  equality  is  found,  the 
position  of  the  Identifying  words  In  the  list  Is  stored  as  this 
Is  the  code  number  of  the  product. 

5-  Construct  Product  Matrices 

Prom  the  list  of  code  numbers  representing  the  gaseous  compounds, 
the  thermodynamic  matrix  Lr  and  part  of  the  A  matrix  are  con¬ 
structed.  Prom  the  list  of  code  numbers  representing  the  con¬ 
densed  compounds,  the  thermodynamic  matrix  Lu*  and  the  remainder 
of  the  A  matrix  are  constructed.  A  melting  point  vector  T^ 
and  a  list  of  the  molecular  weights  of  the  products  are  also 
constructed. 

6.  Calculate  Percent  Composition 

In  this  section  the  Input  composition,  which  can  be  In  either 
moles  or  grams.  Is  normalized  to  100  g.  of  propellant  by  the 
following  method. 

Define  p  =  — - 

E  CiM.  (68) 

1=1 


where  M  =  100  g. 
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Then  the  moles  of  the  1th  reactant  in  the  system  is; 


(69) 


If  is  the  molar  heat  of  formation  of  the  1th  reactant  and 

oJ  Jhe“lac™Sfsyste™'l3  =  °"  “‘^alpy 


P 

r 

1=1 


^1  ^f. 


qJ^nSlJy  ?a“ai™iIJaS"byT 


P  P 

2_  i\  AH^  =  p  f  c  AH. 

iti  *1  A  ^  <^1 


1=1 


(70) 


Se''reactSiL^^?<-®f^®''^®  of  chanps  in  the  heat  of  formation  of 

convenient  to  Include  two  fictitious  reac¬ 
tants  to  the  list  of  propellants.  These  fictitious  reactant-*? 
ave  heats  of  formation  of  +100  kcal./mole  and  -100  kcal  /mole 
ta^ts°to^®®‘  Hence,  by  including  x  moles  of  one  of  these  reac- 

(S?  deo?easSf‘'brioOx’‘JoJ?’  ""S  formation  is  increased 

vor  aecreased;  by  lOOx  kcal.  and  the  mass  remains  unchanged. 


7-  store  Preliminary  Data 


^e  percent  composition,  formulas  and  names  of  the  reactants 
the  formulas  of  all  the  products  as  well  as  Identify*  rSn 

DOrtlon^of  the^’c'*  t®P®-  thereby  allowing  ^large 

portion  of  the  computer  memory  to  be  utilized  again.  • 

8. 


Calculate  a  Starting  Composition 


Semen?  tl^lTenXT'°^  <=he  1th 


P 

T. 


a.  r . 
1=1  1 


(71) 


<V. 


n 
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o 


Inhere  aij  Is 
1th  reactant. 


the  number  of  atoms  of  element  j  per  molecule  -of 
Substitution  of  equation  (69)  Into  (7I)  gives; 


P 

pI 


1=1 


(71) 


number  of  atoms  of  element  <  per  molecule  of 
Ith-reacllant.  Substitution  of  equation  (69)’^iSto  (71)  gives : 


P 

=  p  E  HjH 

1=1 


(72) 


or  In  matrix  notation: 


B  =  pAp  c 


(73) 

1 . e . ,  the 
calculated 


Example  : 

For  the  reaction  mixture; 


NaH4  +  2Pa  - >  products. 


1 

1  corresponds  to  N2H4 

2  corresponds  to  Pa 


1 

2 

3 


1 

corresponds  to  H 
corresponds  to  N 
corresponds  to  P 


The  formula  numbers 
tune,  Cj^j  and  the  mo 
given  below. 


,  the  coefficients  of  the  reaction  mix- 

lecular  weights  of  the  reactants,  Mj^,  are 


H  N  P _ 

1  Cl  Ml 

N2H4  1  4  2  0  1  32.048 

^2  2  0  0  2  2  38.00 

Prom  equation  (68). 

n  ^  _ _ 100.0 

P  32'.o48  +  2  (58rooT 
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Equation  (75)  yields: 


r  1 

4o' 

fl  ] 

4 

[bJ  =  p 

20 

Lo2_ 

[sj  =  P 

2 

4 

B  m 

that  is,  the  system  contains  4p  gram-atoms  of  H,  2p  gram-atoms 
of  N,  and  4p  gram-atoms  of  P.  The  total  v/eight  of  the  system 
is  100  g. 

The  White,  Johnson,  Dantzig  Method  requires  a  starting  composition 
vector  Y  v/hich  is  compatible  v/ith  the  element  balance.  Let  Y® 
be  an  arbitrary  Y  v/hich  will  be  corrected  by  Y^  so  as  to  make 
the  resulting  Y*  compatible  with  the  element  balance; 

L  aijyi*  =  bj  (74) 

1=1 

Substitution  of: 

Yi*  =  y®^  +  y^c  (75) 

into  (74)  produces: 

i=l  J=1 

It  is  desired  to  determine  the  yi*^'s  from  Equation  {76).  Since 
there  are  but  m  equations,  only  m  of  the  yj's  can  be  determined 
from  the  above  equations.  Let  values  of  zero,  therefore,  be 
assigned  to  n-m  of  so  that: 

n  m 

E  ^ijyi°  +  E  aijyic  =  bj  (77) 

1=1  i=l 


In  matrix  notation: 

A'  Y°  +  Ag'  yc  =  B  (78) 

where  As'  is  a  square  m  x  m  matrix  whose  columns  correspond  to 
m  "key"  compounds. 

Solving  (78).  for  yields; 

^  YC  =  [Ag']“^  [B  -^A'Y°]  (79) 
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where  the  elements  of  correspond  to  the  m  k<5y  compounds. 

From  (75)  and  (79)  we  obtain; 

Y  =  Y“  +  [As']‘^  [B  -  A*Y“]  (80) 

v/here  it  is  understood  that  the  corrections  are  added  only  to 
the  y^  's  of  the  key  compounds.  Hence,  given  an  arbitrary  com¬ 
position  vector  Y  and  m  key  compounds  we  can  calculate  a  com¬ 
position  vector  Y*  which  is  compatible  with  the  element  balance 
equations . 

A  convenient  choice  of  Y®  is  a  vector  each  element  of  v/hlch  is 
0.0001  moles.  A  convenient  choice  of  key  compounds  are  the 
chemical  elements.  This  choice  of  "key  compounds”  eliminates 
the  possibility  of  negative  yi*'s. 

Vmite,  Johnson,  and  Dantzig  have  stated  that  a  suitable  starting 
composition.  Y®,  must  satisfy  the  mass  balance  constraints,  l.e.. 
Equation  (7).  They  also  state  that  round-off  errors  in  the 
composition  Introduced  during  the  course  of  the  calculations  can 
be  reduced  by  using  an  equation  similar  to  (8).  These  remarks 
suggest  that  compositions  v/hich  are  not  compatible  v/ith  the 
mass  balance  constraints  have  a  tendency  to  make  the  calculated 
solutions,  l.e.,  X's,  drift  from  the  true  solution. 

It  has  been  found,  however,  that  the  convergence  method,  at  least 
for  a  computer  with  floating  point  hardware,  is  self -stabilizing. 
In  other  words,  the  method  will  seek  out  the  true  solution 
v/hether  the  input  composition  satisfies  the  mass  balance  con¬ 
straints  or  not.  To  verify  this  point  the  following  numerical 
experiment  was  performed  on  a  few  problems.  First,  the  equili¬ 
brium  composition  was  calculated  from  a  starting  composition 
which  agreed  with  the  mass  balance  constraints.  Presumably 
this  equilibrium  composition  was  correct.  The  same  calculation 
was  then  repeated  except  that  the  starting  composition  disagreed 
drastically  with  the  mass  balance  equations.  The  equilibrium 
composition  calculated  in  the  second  case  agreed  with  the  correct 
equilibrium  composition. 

Therefore,  finding  a  starting  composition  which  is  compatible 
with  the  mass  balance  constraints  seems  to  be  an  unnecessary 
operation,  although  composition  convergence  Is  attained  much 
faster  If  a  "correct"  starting  composition  Is  used. 

9.  Chamber  and  Exhaust  Solutions 


The  methods  used  to  obtain  both  the  chamber  and  exhaust  solutions 
have  been  extensively  discussed  In  Sections  II  and  III.  Figure 
2  represents  a  flow  diagram  of  the  method  used  for  computing 
an  equilibrium  composition. 


10.  storing  of  Results 


All  results  are  stored  on  magnetic  tape  as  they  are  obtained. 
After  the  problem  has  been  completed^  a  separate  program  takes 
the  data  from  the  tape  and  prints  the  results.  Since  the  pro¬ 
gram  is  automatic^  this  method  of  storing  the  solution  saves 
time  by  the  nature  of  the  speed  of  writing  on  magnetic  tape 
versus  punching  out  the  ansv;ers .  A  sample  output  is  included 
in  Appendix  II. 

B.  Capabilities  of  the  Computer  Program 

A  program  for  the  Burroughs  220  to  calculate  the  specific  impulse 
has  been  v;rltten  according  to  the  method  lllustustrated  in  Figure 
5-  The  capabilities  of  this  program  will  now  be  considered. 

A  study  of  Figure  2  reveals  that  the  last  operation  of  a  given 
problem  reads  in  the  succeeding  problem  thereby  making  the 
program  automatic.  Any  number  of  problems  may  be  run  without 
manual  intervention.  Probable  error  stops  and  non-converging 
solutions  are  accounted  for  by  a  printout,  which  identifies  the 
difficulty,  follov;ed  by  an  instruction  to  read  in  the  next  pro¬ 
blem.  This  scheme  enables  unattended  overnight  operation  of 
the  program. 

Provided  that  sufficient  thermodynamic  data  are  availabe,  the 
program  is  general.  Computer  storage  capacity,  however,  imposes 
certain  restrictions  on ‘the  size  of  systems  that  can  be  treated. 
These  restrictions  are: 

(a)  14  or  less  reactants  in  the  fuel  mixture. 

(b)  12  or  less  elements  in  the  fuel. 

(c)  12  or  less  condensed  species  present  in  the  products 
or  reaction. 

(d)  90  or  less  products  of  reaction. 

To  date  all  of  the  problems  attempted  (over  6000)  have  remained 
within  these  restrictions. 

The  computer  time  required  to  solve  one  complete  problem  varies 
from  3  minutes,  for  simpler  problems,  to  about  50  minutes  for 
the  most  difficult  problems.  Most  problems  average  about  25 
minutes  of  computer  time. 

Input  data  for  a  given  problem  are  kept  to  a  minimum.  The 
selection  and  preparation  of  thermodynamic  data,  as  well  as  all 
preliminary  calculations  (such  as  calculating  the  gram-atoms  of 
each  element  in  the  fuel  and  enthalpy  of  the  fuel)  are  performed 
by  the  computer.  The  only  input  data  required  are  the  composi¬ 
tion  of  the  reaction  mixture,  the  pressures  of  the  chamber  and 
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nozzle  exit,  and  a  rough  estimate  of  the  chamber  temperature. 

All  other  pertinent  data  are  retrieved  from  magnetic  tape  where 
they  are  permanently  stored.  The  composition  of  the  reaction 
mixture  may  be  entered  either  in  terms  of  the  fuel  compounds  (which 
is  usually  most  desirable),  or  in  terms  of  the  elements  and  the 
fuel  enthalpy. 

An  example  of  input  and  output  from  the  computer  is  given  in 
Appendix  III. 
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VI.  NOMENCLATURE 


A 


Ac 


Ag 

Ap 


-  an  m  X  n  matrix  whose  elements  are  a±j  for  the  n 
products  of  combustion. 

-  an  m  X  nc  matrix  whose  elements  are  aij  for  the 
n^  condensed  products  of  combustion. 

-  an  m  X  ng  matrix  whose  elements  are  aij  for  the 
ng  gaseous  products  of  combustion. 

-  an  m  X  p  matrix  whose  elements  are  a^j  for  the 
p  propellants. 


As 


an  m  X  m  matrix  whose  elements  are  a^j  for  the 
m  "key"  components  used  in  making  a  mass  balance. 


B  -  an  m  X  1  column  matrix  v;hose  elements  are  bj. 

Cp(T)  -  the  heat  capacity  at  constant  pressure  as  a  function 

of  temperature. 


Cpl 


Pi 

P(T) 

Pr(X) 

y- 

3-c 

Ji 

G 

H(T) 


-  the  heat  capacity  at  constant  pressure  at  the 
temperature  cal./mole/“K. 

-  the  free  energy  of  the  1th  species,  cal. /mole. 

-  the  molar  free  energy  of  a  compound  as  a  function 
of  temperature,  in  cal. 

-  the  free  energy  of  a  system  at  a  given  temperature 
as  a  function  of  the  composition  X  divided  by  RT. 

-  an  n  X  1  column  matrix  whose  elements  are 

-  an  nc  X  1  colvimn  matrix  whose  elements  are  the 
for  the  condensed  species. 

-  the  ideal  partial  molal  free  energy  of  the  1th 
species  divided  by  RT. 

-  an  ng  X  1  coliamn  matrix  whose  elements  are  the  ^ 
for  the  rig  gaseous  species. 

-  a  function  defined  by  Equation  (8). 

-  the  enthalpy  of  the  system  as  a  function  of 
temperature,  in  cal. 


-50- 


H(T,  X) 
H(T,  X,  P) 


He(T“) 


AHfj 


-  the  enthalpy  of  the  system  as  a  function  of 
temperature  and  composition,  in  cal. 

I  -  the  enthalpy  of  the  system  as  a  function  of 
temperature  composition  and  pressure,  in  cal. 

the  enthalpy  of  m  grams  of  products  in  the 
rocket  chamber,  in  cal. 

the  enthalpy  of  »i  grams  of  products  at  the 
rocket  exit,  in  cal. 

the  enthalpy  of  the  elements  at  the  reference 
temperature,  cal. /mole. 

-  the  molar  enthalpy  of  the  ith  species,  in  cal./mole. 

the  enthalpy  of  the  system  calculated  at  the 
temperature,  Tn,  in  cal. 

'  the  reactant  temperature 

of  the  Jth  reactant,  in  cal./mole. 

the  column  matrix  corresponding  to  *?. 

'  temperature  which  is  used  to  describe 

the  enthalpy  of  a  system,  units  depend  on  J. 

the  enthalpy  thermal  vector;  a  vector  whose 
elements  are  Ij. 

the  specific  Impulse,  in  seconds. 

"  of  temperature  which  is  used  to  describe 

the  entropy  of  a  system,  units  depend  on  j. 

the  column  matrix  corresponding  to 

-  the  entropy  thermal  vector;  a  vector  whose  elements 
are  Jj. 

the  column  matrix  corresponding  to  1^. 

-  K/RT 

“  temperature  which  Is  used  to  describe 

the  free  energy  of  a  system,  units  depend  on  J. 
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-  the  free  energy  thermal  vector;  a  vector  whose 
elements  are  Kj. 

-  an  n  X  (Mf2)  matrix  known  as  the  thermodynamic 
property  matrix;  its  rows  are  the  T's  for  all 

the  compounds  in  the  system  at  a  given  temperature. 

-  an  nc  X  (N+2)  matrix;  its  rows  are  the  “if' s  for  the 
nc  condensed  species  at  a  given  temperature. 

-  a  5u  X  (N+2)  matrix;  its  rows  are  all  of  the*^'s 
for  the  u  condensable  compounds. 

-  a  r  X  (J»H-2)  matrix;  its  rows  are  the  "t's  for  the 
always  gaseous  compounds. 

-  an  ng  X  (N+2)  matrix;  its  rov;s  are  then's  for  the 
ng  gaseous  compounds  at  a  given  temperature. 

-  the  molecular  v/eight  of  the  ith  species,  gm./mole. 

-  an  n  X  1  column  matrix  v/hose  elements  are  . 

-  the  niamber  of  tenns  used  to  characterize  the 
functional  dependence  of  Cp  v/ith  temperature. 

-  the  pressure  of  the  system,  in  atm. 

-  the  pressure  at  the  nozzle  exit,  in  atm. 

-  the  pressure  in  the  rocket  chamber,  in  atm. 

-  an  n  X  1  coliimn  matrix  v/hose  elements  are 

-  the  ideal  entropy  of  mixing  of  the  1th  species 
divided  by  R. 

-  an  nc  X  1  column  matrix  v/hose  elements  are  the 
entropy  of  mixing  of  the  ng  condensed  species 
divided  by  R.  This  matrix  is  null. 

-  an  ng  X  1  colvimn  matrix  whose  elements  are  the 
for  ng  gaseous  species. 


“  the  quadratic  approximation  to  the  free  energy 
resulting  from  a  Taylor's  Expansion  of  P(X) . 

R  -  the  gas  constant  per  mole,  cal./nole/®K. 

S(T,X,P)  -  the  entropy  of  the  system  as  a  function  of 

temperature,  composition,  and  pressure,  in  e.u. 

%  ■  the  molar  entropy  of  the  ith  product  species,  in  e.u. 

^1  “  the  partial  molar  enthalpy  of  the  1th  product  species. 

In  e.u. 


Sc 


ASi 

T 

Tn 

Tb 

X 

Xc 

X  ° 
Ac 

Xe 

Xe° 

Xg 

Xnew 


-  the  entropy  of  grams  of  products  In  the  rocket 
chamber,  in  e.u. 

-  the  Ideal  entropy  of  mixing  of  the  ith  product 
species.  In  e.u. 

-  the  temperature,  “K. 

-  the  nth  approximation  to  a  solution  temperature. 

-  a  base  temperature,  ®K. 

-  the  reference  temperature,  “K. 

-  an  n  X  1  vector  whose  components  are  x^. 

-  an  nc  X  1  colvimn  matrix  whose  components  are  the 
xi  for  the  nc  condensed  species. 

-  the  equilibrium  X  at  the  chamber  temperature. 

-  the  X  of  the  product  species  at  the  rocket  exit. 

-  the  equlllbriiun  X  at  the  rocket  exit  temperature 
for  equlllbrivun  flow. 

-  a  ng  X  1  column  matrix  whose  elements  are  the  Xj^ 
for  the  ng  gaseous  species. 

-  y  +  XA 


-  an  n  X  1  column  matrix  whose  components  are  the  yj_. 

-  an  n  X  1  column  matrix  whose  components,  yj^*, 
satisfy  the  mass  balance  constraints. 
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yO 

Yg 

YgD 


-  an  arbitrary  n  x  1  column  matrix  whose  elements 
are  yj^®. 

-  an  n  X  1  column  matrix  whose  components  are  y^®* 

-  an  ng  X  1  column  matrix  whose  elements  are  the  yj^ 
for  Che  ng  gaseous  species. 

-  an  ng  X  ng  diagonal  matrix  whose  diagonal  elements 
are  the  ^y^  for  the  ng  gaseous  species. 


a 


ai 


I 


b 

c 

Ci 

c(T) 

d 

e 

fj(T) 

Ij 

1 


m 


-  an  empirical  constant  used  to  describe  heat  capacity 
data. 


-  the  formula  designation  vector;  a  vector  whose 
elements  are  a^i,  &i2>  • • ‘^im* 

-  a  formula  designation  number  v/hich  indicates  the 
number  of  atoms  of  element  j  in  a  molecule  of 
species  i. 


-  an  empirical  constant. 

-  the  total  nvimber  of  gm.  atomic  weights  of  element 
J  present  in  the  mixture. 

1/2 

-  a  constant;  in  Equation  (l)  c  =  9*5281 

-  the  number  of  moles  of  the  ith  reactant  in  the 
reaction  mixture  (propellant). 


-  the  number  of  condensed  compounds  in  a  given  system 
at  the  temperature  T. 


-  a  constant. 


-  a  constant. 

-  a  funct.'^on  of  temperature  used  to  describe  constant 
pressure  heat  capacity  data. 

-  empirical  constant  used  to  describe  thermodynamic 
properties,  see  Equations  56,  45  and  44. 

-  a  thermodynamic  property  vector;  a  vector  whose 
elements  are  Ij. 

-  the  number  of  chemical  elements  in  a  given  system. 
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r'l 


n 

"c 

n 

P 

r 


g 


u 

v(T) 


X 

y 


Ys 

y 

n 

a 

yi 

yi' 

yi^ 


-  the  total  number  of  product  species  of  a  system. 

-  the  number  of  condensed  product  species  of  the 
system  at  a  given  temperature. 

-  the  n\imber  of  gaseous  product  species  of  the  system  . 
at  a  given  temperature. 

-  the  nvmiber  of  reactants  (propellants). 

-  the  number  of  compounds  in  a  given  system  v/hich  are 
alv;ays  gaseous  above  1000“K. 

-  the  nvmiber  of  moles  of  the  Jth  reactant  in  the 
propellant  mixture. 

-  the  number  of  condensable  compounds  in  a  given  system, 

-  the  nvmiber  of  vaporized  compounds  in  a  given  system 
at  the  temperature  T. 

-  the  svimmation  of  the  gaseous  x^. 

-  an  np.  X  1  column  matrix  v/hose  elements  are  all  x/  . 

®  y 

-  the  corrected  value  of  x  given  by  (26). 

-  the  number  of  moles  of  the  1th  product  species; 
x^  Is  a  variable. 

-  the  number  of  moles  of  the  sth  product  species  where 
s  Identifies  the  smallest  X,  (see  Equation  24). 

-  the  yj^  which  corresponds  to  Xq. 

-  the  summation  of  the  gaseous  y^. 

-  the  number  of  moles  of  the  1th  product  species; 

Is  one  of  the  coordinates  of  Y,  a  fixed  point 
about  which  a  Taylor's  Expansion  Is  made. 

-  the  niimber  of  moles  of  the  1th  species;  an  element 
of  Y*. 

-  an  arbitrary  number  of  moles  of  the  1th  species. 

-  Yi*  -  yi°. 


lit 
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-  an  m  X  1  column  matrix  v/hose  elements  are  Oj. 

-  the  atomic  weight  of  the  Jth  element. 

-  a  scale  factor  which  reduces  the  step  length  given 
by  A.  See  Equation  (  23) . 

-  the  largest  X  which  prevents  x^from  becoming 
negative . 

-  the  smallest  of  all  the 

-  the  mass  of  the  propellant  system,  in  gms. 

-  an  m  X  1  vector  v/hose  elements  are  the 

-  a  Lagrange  multiplier  associated  with  the  mass 
balance  equation  for  the  Jth  element. 

-  a  vector  whose  elements  are  the  itj. 

-  the  number  of  reaction  weights  In  a  system  of 
mass  u. 

-  a  constant  multiplier. 

-  a  factor  used  in  estimating  Improved  temperatures 
when  enthalpy  is  used  in  the  convergence  criterion. 

-  a  factor  used  in  estimating  improved  temperatures 
v/hen  entropy  is  used  in  the  convergence  criterion. 

-  an  n  X  1  column  matrix  or  vector  which  indicates 
the  directions  toward  an  improved  solution  X-Y. 

-  an  element  of  A,  x^  -  y^. 

Subscripts 

-  chamber  conditions;  condensed  species. 

-  exit  conditions . 

-  gaseous  species. 

-  a  dummy  index  usually  used  to  designate  a  species. 

-  a  dummy  index  usually  used  to  designate  a  chemical 
element . 


a  dummy  Index. 

Superscripts 

denotes  the  transpose  of  a  matrix, 
denotes  the  Inverse  of  a  matrix. 
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VIII.  APPENDICES 


APPENDIX  I 

THE  THERTdODYNAMIC  PROPERTY  VECTORS 
OP  THE  COMBUSTION  PRODUCTS - 


The  general  formulation  of  the  thermodynamic  property  vectors 
given  in  Section  IV  is  now  applied  using  the  heat  capacity 
expression  suggested  by  Thomson'*  and  the  thermodynamic  property 
data  compiled  by  The  Dow  Thermal  Laboratory.®  The  thermodynamic 
property  vectors  resulting  from  the  above  source  are  tabulated 
at  the  end  of  this  Appendix. 

The  constant  pressure  heat  capacity  expression  proposed  bv 
Thomson  is  of  the  form: 


Cp  =  li  +  I2  +  I3  +  1,^  T®  +  I5  T-2  (81) 

where  T  is  the  temperature  in  “K. ,  divided  by  1000  and  the  li's 
are  constants  to  be  determined  from  experimental  measurements. 


Follov/ing  the  method  prescribed  by  Thomson,  the  constants  li 
through  I5  v/ere  calculated  from  values  of  the  heat  capacity 
measured  at  five  characteristic  temperatures.  If  Ti  to  T5 
and  Cpi  to  Cps  are  the  characteristic  temperatures  and  their 
corresponding  heat  capacities,  respectively,  then  (8I)  leads  to 
the  following  matrix  equation: 


1  Ti  Tf  Tf  tI^ 

"li 

Cpa 

1  Ta  Ti  Ti  Ti® 

la 

Cp3 

= 

1  T3  Tf  Ti  Ti® 

I3 

Cp4 

1  T4  Tf  Tf  Tf® 

I4 

1  Ts  Tf  Ti  Ti® 

,15 

(82) 


Solving  for  unknov/ns  lj_'s  gives; 


li’ 

-1 

Cpi 

la 

I3 

crj" 

Cpa 

Cp3 

I4 

Cp4 

fp3 

(85) 


where  l  is  the  square  matrix  of  Equation  (82).  For  species  not 
undergoing  phase  changes,  the  five  characteristic  temperatures 
are : 


1000,  1600,  3000,  4500,  and  5700°K. 


Howeverj  If  a  phase  change  did  occur,  then  five  temperatures 
with  spaclngs  proportional  to  those  above  were  chosen  so  as  to 
cover  the  complete  temperature  range  of  each  phase.  (Table  III 
contains  the  temperature  range  used  In  calculating  the  li's 
for  each  compound.)  Equation  (83),  therefore.  Is  used  to  cal¬ 
culate  li  through  I5. 

The  quantities  la  and  I7  are  calculated  according  to  (43a)  and 
(44a): 

5 

.  la  =  H(Tb)  -  J  Ij  Ij(Tb)  (43a) 

J=1 

I7  =  S(Tb)  -  f  Ij  Jj(Tb)  (44a) 

J=1 

Before  these  expressions  can  be  applied,  the  functions  H,  S,  I, 
and  J  must  be  evaluated  at  T  =  Tb- 

Absolute  entropies  are  generally  obtainable,  but  absolute 
enthalpies  are  not.*  Consequently,  S(Tb),  In  (44a),  Is  the 
absolute  entropy  at  the  base  temperature,  Tb,  while  H(Tb),  In 
(43a),  must  be  a  relative  enthalpy.  A  convenient  reference 
state  for  evaluating  enthalpies  Is  the  elements  at  standard 
conditions.  If  this  reference  state  Is  chosen  then: 

H(Tb)  =  (H(Tb)  -  H(T*)]  +  (H(T“;)  -  He(T“))  (74) 

Both  of  the  bracketed,  terms  on  the  right  side  of  (74)  are 
available  from  the  Thermodynamic  Tables®.  The  first  term  Is 
the  change  In  enthalpy  In  going  from  the  reference  temperature, 
T®,  to  the  base  temperature,  Tp,  while  the  second  term  Is  the 
heat  of  formation  from  the  elements  at  T®.  The  entropy  at  the 
base  temperature,  S(Tb),  Is  also  available  from  the  Thermodynamic 
Tables. ® 

The  functions  I  and  J  are  defined  by  Equations  (39)  and  (40). 


Ij(T)  =  ^j(T)  dT  (39) 


The  above  Integrations  are  straightforward  and  are  given  In 
Table  II. 

* 

This  stems  from  difficulties  In  measuring  the  zero  point  energy. 

^  This  reference  state  Is  especially  useful  In  treating  multi- 
component  systems  since  It  Is  Independent  of  the  composition. 
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TABLE  II 


Elements  of  the  Thernial  Vec tors 


J 

fj(T) 

Ij(T) 

Jj(T) 

Kj(T) 

1 

1 

T 

In  T 

T(l-ln  T) 

2 

T 

1/2  T® 

T 

-  1/2  T® 

3 

T2 

1/5  T® 

1/2  T^ 

-  1/6  T® 

4 

qi3 

1/4  T'* 

1/5  T® 

-  1/12  T* 

5 

T-2 

-T-i 

-  1/2  T‘2 

-  5/2  T"^ 

6 

- 

1 

0 

1 

7 

- 

0 

1 

-T 

Hence,  from  H(Tb)  calculated  via  (8if),  S(Tb)  obtained  from  the 
Thermod^amlc  Tables®,  and  Ij(Tb)  and  JjiTb)  calculated  from 
Table  II,  Iq  and  I7  can  be  evaluated  using  (45a)  and  (44a). 


of  Equation  (51)  or  (52)  reveals  that  Tb  must  be  chosen 
within  cne  temperature  range  for  which  li  through  I5  are  valid. 
Choosing  a  Tb  outside  this  range  is  satisfactory  only  if  l, 
through  I5  represent  the  behavior  of  Cn  versus  T  to  the  tempera¬ 
ture  Tb«  This  condition,  hov/ever,  does  not  generally  prevail 
^cause  of  slight  variations  between  the  observed  heat  capacity 
data  and  those  predicted  from  li  through  I5,  slight  variations 
in  Is  and  I7  occur  if  T^  is  allov/ed  to  vary.  Therefore,  rather 
than  choose  Is  and  I7  corresponding  to  a  given  Tb,  it  seems 
best  to  calculate  the  average  values  of  Is  and  I7  based  on 

Values  of  Is  and  I7  entered  in  Table  III  are  the 
arithmetic  means  of  Is  and  I7  calculated  at  Tb's  equal  to  the 
five  characteristic  temperatures  described  above. 


By  itself  the  thermodynamic  property  vector  1  has  little  signifi¬ 
cance;  it  must  be  associated  with  the  enthalpy  thermal  vector 
I,  the  entropy  thermal  vector,  j,  or  the  free  energy  thermal  * 
vector,  K,  before  it  takes  on  real  significance.  The  elements 
of  these  thermal  vectors  are  listed  in  Table  II.  Notice  that 
fi  through  fg  are  defined  by  (8I),  and  for  J=6  and  7  are 
defined  by  (45)  and  (44)  and  the  Kj's  are  defined  by  (48). 


Using  the  procedure  outlined  above,  the  elements  (li  through 
j  thermodynamic  property  vectors  of  98  combustion  pro¬ 

ducts  have  been  calculated.  The  results  of  these  calculations 
are  entered  in  Table  III  below. 


Values  of  Cg,  H,  and 
vectors  in  Table  III 
than  1  per  cent . 


S  generated  from  the  thermodynamic  property 
agree  with  the  tabulated  values®  to  better 
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APPK'iDIX  II 


INPUT 

Input  data  for  a  given  problem  consist  of: 

(l)  A  test  word  v^hlch  Indicates  whether  the  input 
is  in  moles  or  grams,  and  whether  a  modifica¬ 
tion  in  the  number  of  product  species  is  to 
be  considered. 

(11 )  The  propellant  composition  designated  by  two 
numbers.  The  first  number  is  the  number  of 
moles  or  grams  of  the  given  reactant  in  the 
reaction  mixture  and  the  second  number 
identifies  the  reactant. 

(ill)  The  chamber  and  exit  pressures. 

(Iv)  A  rough  estimate  of  the  chamber  temperature. 

(v)  An  identifying  run  number. 

(vl)  Any  appropriate  changes  in  the  melting  points 
of  the  condensable  species. 

(vli)  Any  modification  of  the  product  species. 

Figure  5  is  a  sample  input  sheet  for  a  lithium  aluminum  hydride- 
nitronlum  perchlorate-  polyethylene  system.  The  pressures  and 
temperature  are  given  in  floating  decimal  number. 

There  is  also  Included  the  output  data  sheet  obtained  from  the 
input  (Figure  4). 


Pig.  5  -  Sample  Input  for  Specific  Impulse  Calculation 


Run  Number  5736 

No.  Products  58  No.  Gases  5^*  No.  Condensed  Phases 


Percent 

Mole  Weight  Formula 
29. Ho  18.00  LiAlH4-s 
26.^12  62.00  ClNOo-3 
^J^.19  20.00  CaH4-s 


Name 

Lithium  Aluminum  Hydride 
Nltronlum  Perchlorate 
Polyethylene 


APags 


Species 

Chamber 

Moles 

Exhaust 

Moles 

Exhaust 

Grams 

A1 

.001737 

AlCl 

.010421 

AICIO 

.000019 

AlCla 

.003354 

. 000001 

AICI3 

.000008 

AlH 

.000458 

AlHO 

.000021 

AlHOa 

. 001975 

AIN 

.000001 

AlO 

.000655 

AlaO 

.002284 

AlaOa 

.000022 

HCN 

.000156 

. 000005 

. 000086 

CHNO 

.000007 

.000005 

CHO 

.001285 

. 000009 

. 000264 

CHaO 

.000021 

. 000012 

CH4 

.000002 

. 000007 

ON 

. 000006 

CO 

1.594454 

1.570447 

58.587588 

COa 

. 029877 

.055350 

2.455122 

Cl 

. 005279 

. 000005 

. 000101 

HCl 

.075704 

.008211 

.299418 

HOCl 

. 000001 

LiCl 

. 528966 

. 402002 

17.045689 

Cla 

. 000001 

LlaCla 

. 000544 

.007995 

. 677887 

LI3CI3 

. 000002 

. 600208 

H 

. 125172 

. 000919 

. 000926 

LIH 

.009217 

. 000107 

. 000850 

LiOH 

. 082462 

.054228 

1.298659 

Pig.  4  -  Output  Data 
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Sheet 


Run  Number  5756  (Continued) 


Species 


NH 

OH 

Ha 

HaO 

LI 

LIN 

LIO 

LlaO 

N 

NO 

Na 

0 

Oa 

AlaOa-s , 1 


Chamber 

Moles 


. 000057 
.008881 
1.897070 
.525619 
.0^886^1 
.0050^1 
.000'I22 

.000108 

.000008 

.000361 

.211302 

.000335 

.000025 

.225550 


Exhaust 

Moles 


.000009 

2.052152 

.510^85 

.001925 

.000016 

.000002 


.215100 


.2571^«2 


Exhaust 

Grams 


.000156 

096817 

5.595658 

.015561 

.000555 

.ooooon 

.000064 

.000002 

5.970205 


24 . 178961 


Temperature,  °K. 
Pressure,  atmos. 
Enthalpy,  kcal. 
Entropy,  e.u. 
Moles  of  gas 
Avg.  mol.  wt.,  g. 
Gamma 

Velocity,  cm./sec 
Impulse,  sec. 


3278.838  1781.701 

68.046  1.000 

-17.288  -113.398 

267.536  267.542 

4.566  4.457 

16.864  17.011 

1.271  1.297 

283601.740 
289 . 182 


Pig.  4  (Contd.) 
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